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Abstract. Today, the electron affinity is experimentally well known for most 
of the elements and is a useful guideline for developing ab initio computational 
methods. However, the measurements of isotope shifts on the electron affinity 
are limited by both resolution and sensitivity. In this context, theory eventually 
contributes to the knowledge and understanding of atomic structures, even though 
correlation plays a dominant role in negative ions properties and, particularly, in 
the calculation of the specific mass shift contribution. The present study solves the 
longstanding discrepancy between calculated and measured specific mass shifts on 
the electron affinity of chlorine [Phys. Rev. A 51, 231 (1995)]. 
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1. Introduction 

Despite the difficulty to measure electron-affinities ( e A), experimental values can 
be considered as "exact" by theoreticians for most of the elements. Due to their 
remarkable properties as weakly bound atomic systems largely governed by electron 
correlation, negative ions have always attracted broad attention from chemists and 
atomic physicists pQ. With this respect, the field of electron affinities constitutes a 
playground of choice for testing the computational methods [2] [3] . The estimation of 
the mass polarization operator expectation value that partly makes the level isotope 
shift (IS), also highly sensitive to correlation 0], offers other stringent tests for 
assessing the reliability of theoretical methods. Unfortunately, the mass shifts of 
transition frequencies or electron detachment thresholds are rather difficult to extract 
from experimental studies. Limitations arise from the size of isotope shifts and/or 
from the low abundances of the isotopes to be studied, in which cases isotopically 
enriched samples [5] or in-source spectroscopy [Bj could be used. One way of deducing 
a pure mass shift value from observation is the use of King plots [7] |S], requiring 
measurements of isotope shifts on a transition for at least three pairs of isotopes. The 
nuclear stability valley narrowing for low atomic numbers, there are less possibilities 
of building the required experimental King plots when approaching the domain of 
systems for which it is realistic to do highly correlated calculations. Any other way 
to estimate the field and mass shift contributions separately must contain theoretical 
inputs [91 [TO]. In cases where either the field shift (FS) or the mass shift (MS) is a 
small contribution, the error on the theoretical part is often neglected, although it 
involves in some cases an unknown systematical error in the data analysis . 

In the context of the physics of highly-charged ions, the balance between the mass 
and field shift contributions to the level and transition isotope shifts is worthwhile to 
investigate systematically along an isoelectronic sequence [TO] for which the increasing 
importance of a proper relativistic treatment towards the high-Z region becomes 
obvious P31 HU [15]. However, for the study of the isotope shifts on the electronic 
structure of negative and neutral atoms of the second and third periods of the periodic 
table, electron correlation still largely dominates the relativistic corrections and non- 
relativistic approaches constitute a good starting point to investigate what makes an 
isotope shift "normal" or "anomalous" [TOJ [T7] . 

In the field of negative ions, the laser photodetachment threshold technique was 
developed and successfully applied to the measurement of electron affinities [18] and 
of the detachment thresholds corresponding to different fine-structure levels of the 
negative ion and the neutral atom jTSJ [20] . In this method, the onset of the photode- 
tachment process is measured directly using a tunable laser source. In combination 
with improved methods for the measurement of the photoelectron kinetic energies such 
as photodetachment microscopy [21 and slow-electron velocity-map imaging [22], the 
technique of laser photodetachment at threshold represents the most precise way of 
determining photodetachment thresholds experimentally [23) . Amongst the recent ap- 
plications of photodetachment microscopy, let us cite the first experiment realized in 
phosphorus |24j with the excitation of the parent neutral atom out of the fundamental 
spectral term, and the measurement of the electron affinity of selenium with an ac- 
curacy of 1 fieV [25] . The possibility of applying the tunable laser photodetachment 
spectroscopy to the measurement of the isotope shift on the electron affinity has been 
demonstrated for the first time by Bcrzinsh et al J35] for the 35 > 37 C1 isotopes. The 
photodetachment microscopy technique has been applied later to measure the electron 
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affinities of 16 and 18 O separately from a natural sample [27] and deduce the IS on 
the oxygen electron affinities for 16,18 0. This work has been extended to 17 [25] with 
the assistance of theoretical calculations for the estimation of the hyperfme structures. 
The sulfur electron affinities were measured more recently by photodetachment mi- 
croscopy for the two isotopes 32 S and 34 S demonstrating the ability of ab initio 
methods for estimating the isotope shift on the electron affinity. 

More than fifteen years ago, measurements and many-body calculations were 
reported for the isotope shift in the chlorine electron affinity in the pioneer work 
of Bcrzinsh et al [26] . The theory-experiment agreement found for the specific 
mass shift (SMS) was satisfactory as far as the order of magnitude is concerned 
but a serious discrepancy lay in its sign. Still, the obtained resolution by tunable- 
laser photodetachment spectroscopy was remarkable. More recently, the isotope 
shift on the electron affinity of sulfur was investigated, both experimentally and 
theoretically It was established that large scale closed-core multiconfiguration 

Hartree-Fock (MCHF) calculations lead to a reasonable theory-experiment agreement 
but it was also shown that core effects are large and cannot be neglected. The 
success found in 29J for the S/S~ systems incited the present authors to revisit 
theoretically the IS on the electron affinity ( e A) of chlorine, involving only one more 
electron in the neutral/negative ion balance. Section [2] presents the needed theoretical 
background, introducing the mass and field isotope shift on the electron affinity. 
The correlation models and construction of the configuration spaces are described in 
section [3j Section |4~1~1 and section |4~2"] present, respectively, the results of the valence 
correlation models and of the open-core configuration interaction (CI) approach. 
The final theoretical results are summarized and compared with Berzinsh et al 's 
experimental and many-body calculations results [2l)j in section [5] 

2. Electron-affinity and isotope shifts 

2.1. The experimental electron affinity as a guideline 

Like in our previous work on the isotope shift in the sulfur electron affinity [29] . the 
experimental electron affinity 

% exp = E(Cl 3s 2 3 P 5 2 P° /2 ) - E(Cr 3s 2 3 P 6 1 S ) (1) 

is used as a precious guideline to set efficiently pathways in the variational 
configuration spaces. Both 35 C1 and 37 C1 isotopes have a nuclear spin / = 3/2. 
For both isotopes, the ground state of the neutral chlorine is therefore split into four 
hyperfme structure levels F = 0, 1, 2, 3 where the F — level has the lowest energy, 
while the ground state of the negative ion does not show hyperfme structure since 
J = 0. The electron affinity measured by Berzinh et al [26) [Tj. 

A oxp = 29 138.59(22) cm" 1 = 3.612 724(27) eV 

is defined from the threshold energy and thus strictly corresponds to the difference 
between the lowest hyperfme level of the ground state of the neutral and the ground 
state of the anion. This represents a difference with eq. ([!]) of only 700 and 586 MHz for 
35 C1 and 37 C1, respectively and is one order of magnitude smaller than the uncertainty 
of the measured electron affinity. In the present work, we adopt a non-relativistic 
variational approach for targeting electron correlation (see Section [3]) and we need 
to estimate the average experimental electron affinity that would be measured if not 
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resolving the fine-structure thresholds due to the J-splitting of the chlorine ground 

1/2-3/2- 



term CI 3s 2 3p 5 2 P1i2-MT This average electron affinity can be expressed as 



e A^L = E(C\ 2 P°) - E(cr 1 S) 



= (4E 3/2 + 2E 1/2 ) _ 
6 

_e A , gl/2 ~ E 3/2 

Using the NIST fine structure energy separation (E 1 / 2 — -E3/2 = +882.35f5 cm -1 ), 
one finds 54^ = 29 432.70(22) cm -1 = 3.649 189 cV, from which one substracts the 
non-fine structure contribution Ai? NF = —0.01509 eV reported by de Oliveira et al [2] 
to estimate the reference non-relativistic electron affinity 

- AE NF = 3.664 28 eV = 0.134 660 E h . (3) 

Adopting the (A' > A) convention where A is the mass number, the isotope shift 
(IS) on the A is defined as 

IS(A', A) = SA = e A(A') - e A{A) . (4) 

It can be expressed as the sum of the mass shift (MS) and field shift (FS) contributions 

6 e A = S e A MS + 5 e A FS . (5) 

2.2. Isotope Mass Shift 

The mass shift of an atomic energy level is the energy displacement due to the inclusion 
of the dynamics of the nucleus in the Hamiltonian. The proper non-relativistic 
quantum mechanics treatment of the separation of the center of mass motion and 
the motion of the N electrons relative to the nucleus can be found in Johnson [30] and 
Bransden and Joachain [3T]. A first (one-body) effect of the finite nuclear mass is to 
scale the infinite mass energies by the ratio fi/m e = Ma/ '(Ma + m e ). The shift of the 
energy from the infinite-mass value defines the normal mass shift (NMS) 

= Em — £00 = —-rz — r = — ~VT~ = 1i~f~Em -(6) 

M A + m e M A M a 

A second (two-body) correction to the energy, referred to as the specific mass shift 
(SMS) is given by 

,F MA /f \ 1 /f \ m 

SEs ™ = 2(M A + me y \X> ^ ) = (Ma-TW) m e \g Pi • P 7 • (7) 

Combining ([6| and ^ for estimating the level mass isotopic shift (6 Ems — ^nms + 
3Esms)i h is easy to obtain the expression of the mass contribution to the isotope 
shift on the electron affinity Q 

SAms = <^NMS + <5MgMS (8) 



A5sms 



M< M I \M' + m R M - 

where A5*sms is the difference 

ASsms = S'sms(A) - Ssms(X~) (9) 
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of the expectation values defining the SMS parameters 



SSMS 



N 








i<j 





calculated for the ground states of the neutral atom X and the negative ion X~ 
that j8]| can be shown to be strictly equivalent to 

h 2 



(10) 



Note 



S e A M s = 



Al 



M' 



"A — 



M' 



rn,. 



-AS'sms 



A5sms 



(11) 



where the first term alone corresponds to the prescription of King 7jJ] 

The expectation value ( [l0| is dominated by its Hartree-Fock value i.e., in CI, by 
the (s,p) Vinti integrals [33] which contribute negatively. The differential SMS Q 
between the neutral atom and the negative ion is itself dominated by the Vinti integrals 
involving the additional 3p electron in Cl _ . In general, the mass polarization difference 
ASsms is therefore positive, as expected. The first term appearing in the mass 



shift on the electron affinity (11) explicits the antagonism between the two positive 
contributions, ie. the electron affinity 54oo by itself and the SMS parameter variation 
(H 2 /m e )ASsMSj from which strong cancellation between the electronic part of the 
NMS and SMS can be expected : classically, electrons tend to move in opposite 
directions, as beautifully illustrated by Krause et al U in the two-electron case and 
the nucleus to be slower that expected from the one-electron contributions. The SMS 
term is usually the bottleneck in the calculation of the isotope shift of light systems and 
most of our effort is aimed at a robust estimation of this effect in the non-relativistic 
level of approximation. 

The nuclear masses have been estimated by substracting the total electrons 
mass from the atomic masses of 35 C1 and 37 C1 taken from [51] . ie. M( 35 C1) = 
34.95952682(4) u and M( 37 C1) = 36.95657673(5) u. 



2.3. Isotope Field Shift 

According to 35 36 , the field shift on the electron affinity can be estimated from the 
following expression 



S e A FS = 4-7T 



hn 3 

na ft y\AA 

Jz nz > 



(12) 



where pnr(O) is the spin-less total electron density [37] at the origin, and where f(Z) 
corrects for the fact that we use the non-relativistic electronic density for a point 
nucleus. Its value f(Z = 17) 35 ' 37 = 16.940 mK/fm 2 = 507.8 MHz/fm 2 is taken from 
Aufmuth et al 's compilation [55] , 

The knowledge on the nuclear charge distributions of chlorine isotopes is 
unfortunately limited [10] , The root mean square nuclear charge radii (r 2 ) 1 / 2 resulting 
from elastic electron scattering measurements have been reported by Briscoe et al |39j 
for both isotopes 

1/2 



3.384(15) fm , 
v /35 - 3.388(15) fm , 

X The second term arises from the all-order treatment of the NMS (see also footnote (5) of |32p 



( r /37 
/„2\l/2 
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on the basis of phase-shift fits to a three-parameter Fermi distribution. For 35 C1, 
the latter value strongly disagrees with the {r 2 )\^ — 3.335(18) fm rms-value 
extracted from earlier muonic X-ray measu rem ents |40j . To evaluate the difference 
(5^7 ,2 ) 35,3T = (r 2 ) 37 — ( r2 )35 appearing in |l2|, we choose for consistency, the two 
above values recommended by Briscoe et al "3D and also adopted in De Vries et aFs 
compilation ie. 

S ( r 2 ) 35 ' 37 = -0.03(24) fm 2 . (13) 

Note that Berzinh et al used for their FS estimation [26] the value of 6(r 2 ) 35 ' 37 = 
+0.12(12) fm 2 that can also be derived from Briscoe et al 's charge distribution 
parameters, but obtained using two-parameter Fermi models. This positive value 
is in line with the result S (r 2 > ' = +0.127 fm 2 that can be calculated from the rms 
of Angeli's compilation |42j that have been determined by averaging the charge radii 
extracted from different sources |43j . In view of the large dispersion of the nuclear 
data, we think that it is preferable to use the same reference for both isotopes to 
maintain some coherence. 



3. Computational strategy 

3.1. MCHF calculations 

We use multiconfigurational Hartree-Fock (MCHF) wave functions built on 
orthonormal numerical radial orbitals using the ATSp2k package [U]. Such an 
expansion is written 

* MCHF (rLS) = Ci HuLS) , (14) 

i 

where the configuration state functions (CSF) (&i("fiLS) are symmetry-adapted linear 
combinations of Slater determinants. In MCHF theory, both one-electron radial 
functions {P n i( r )} intervening in the construction of the CSF's and the interaction 
coefficients {c.j} are optimized [151 HBJ, while the orbitals are fixed in configuration 
interaction (CI) calculations [47] , 

3.2. Construction of the CSFs spaces 

First, a multireference configuration set (MR) is selected. Then, the "multireference- 
interacting CSF space" (MR-I) is generated as the set of all CSF's interacting with at 
least one CSF of the MR, ie. 

$ 2 (7 4 LS) E MR-I & 3$ fc e MR with ( 7i LS \H\ lk LS) ^ (15) 

whatever the one-electron radial functions are. For a given n max l max defining the 
orbital active space \n max l max ] of subshells with n < n max and I < l ma x, the MR- 
^\n max l max \ set is a subset of the single and double excitations of the MR (MR- 
SD\n max l max ]). As advocated in Rcf. [22], we choose to couple the subshells in 
decreasing order of n and I. 
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3.3. Valence CSF spaces 

Keeping the [Ne] core (ls 2 2s 2 2p 6 ) closed, we use the following multireference 

MRD = [Nc]{3s, 3p} 6 {3, 4} 2 (16) 

for CT~, allowing two electrons to "leave" the spectroscopic subshells 3s and 3p and 
to be excited in the orbitals of the n = 3,4 shells. Since this list is limited to 
double excitations of the main configuration, the MRD-I expansions contain at most 
quadruple excitations. Heeding the strong 3p 2 — > 3d 2 ,Ap 2 interactions, we also explore 



the MR constructed by merging the MRD ( 16 ) with the single excitations of the two 
configurations 3s 2 3p 4 3d 2 and 3s 2 3p 4 4p 2 in the shells n — 3, 4. This list, denoted 
MRT, allows quintuple excitations in the final MRT-I lists. 

For CI, we adopt two different approaches, leading to significantly different 
{Pni{r)} sets. Firstly, we select only the most important configurations in the 
multireference, omitting the single excitations conserving I 

MR1 4 = [Ne]{3s 2 3p 5 ,3s 2 3p 3 3d 2 ,3s 2 3p 3 4p 2 ,3s 1 3p 5 3d 1 }, (17) 

MR1 5 = MRI4 U {Ss^SdHf 1 } ■ (18) 

Secondly, we use a systematic approach based on the following MR space 

MR2 = [Ne]{3s, 3p} 5 {3} 1 {3, 4} 1 . (19) 

MCHF calculations are performed for various n max l max , ranging from 4/ to Ilk, the 
core orbitals being either optimized (MCHF) or fixed to their HF shape (FC-MCHF). 
Table [T] presents the configuration weights obtained in the FC-MCHF calculations. 
These weigths are defined as the summed contribution of all CSFs <&i belonging to a 
given configuration: 



w 



I E c h (2°) 



Contrary to the MR1 4 (C1), MR1 5 (C1) and MRD (CI - ) lists, the use of MR2(C1) 
and MRT(C1~) reference lists favors the apparition of large weights for the single 
substitutions 3s — > 4s and 3p — > Ap by a mechanism identified by the authors [48 . 
The importance of these configurations comes from the strong interaction of the total 
wave function with the Brillouin state of the main configuration, permitted by a large 
deviation of the orbitals from their HF shapes. Very similar wave functions can be 
constructed with and without predominance of the nl — > n'l single excitations, the 
limit case occurring when the wave function is invariant under the rotations between 
the orbitals of a given I. When these single excitations dominate the expansion, the odd 
3^ — > 41 (I = s,p) excitations - i.e. the substitution of an odd number of 31 electrons by 
4/ electrons - are favored while the even ones are suppressed. In particular, this means 
that the MR2-I model can mimic, thanks to the orbital optimization, the presence of 
the 3p 2 — > 4p 2 excitation in the reference. This comes with the price of relatively 
"dilute" eigenvectors and a slower numerical convergence of the MCHF calculations. 
Other manifestations of this effect can be found in Ref. . 



3.4- Core-Valence CSF spaces 

Following the hierarchy of Table [T] subsets of p configurations, denoted MR p , 
are selected. For each p, all the single and double subshell substitutions of the 
corresponding MR p in [life] are generated, allowing at most one hole in the n = 2 shell 
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Table 1. Weights of the configurations composing the MRf4, MRI5 and MR2 
sets of CI and the MRD and MRT sets of CI - in the corresponding MR-Ifllfc] 
wave functions with the core-orbitals fixed at the HF approximation. # is the 
index of each configuration. 





CI 


Sp 5 2 P° 






MRf4 




MR2 






config. 


w 


config. 


w 


# 


3s 2 3p 5 


0.9569 


3s 2 3p 5 


0.9374 


1 


3s 2 3p 3 3cZ 2 


0.1903 


3s 2 3p 3 3d 2 


0.1835 


2 


3s 1 3p 5 3d 1 


0.1143 


3s 2 3p 4 4p 1 


0.1792 


3 


3s 2 3p 3 4p 2 


0.0831 


3s 1 3p 5 3d 1 


0.1098 


4 






os op *±s 


u.uy iu 


O 






3s 1 3p 4 3d 1 4f 1 


0.0690 


6 


MRf 5 




3s 1 3p 4 3d 1 4p 1 


0.0629 


7 


config. 


w 


3s 2 3p 4 4/ 1 


0.0485 


8 


3s 2 3p b 


0.9562 


3p 5 3d 2 


0.0427 


9 


3s 2 3p 3 3cZ 2 


0.1900 


3s 2 3p 3 3d 1 4(2* 


0.0382 


10 


3s 1 3p 5 3d 1 


0.1140 


3p 6 4p 1 


0.0170 


11 


3s 2 3p 3 4p 2 


0.0831 


3s 1 3p 5 4rf 1 


0.0150 


12 


3s 1 3p 4 3d 1 4/ 1 


0.0704 


3p 5 3ci 1 4s 1 


0.0074 


13 






3p 5 3d 1 4d 1 


0.0062 


14 






3s 2 3p i 3d 1 4s 1 


0.0014 


15 




ci- 


3p 6 1 S 






MRD 




MRT 






config. 


w 


config. 


w 


tF 


3s 2 3p 6 


0.9469 


3s 2 3p 6 


0.9124 


1 


3s 2 3p 4 3cZ 2 


0.2090 


3s 2 3p 5 4p 1 


0.2592 


2 


3s 2 3p 4 4p 2 


0.1283 


3s 2 3p 4 3d 2 


0.2066 


3 


3s 1 3p 5 4s 1 4p 1 


0.0737 


3s 2 3p 4 4p 2 


0.0876 


4 


3s 1 3p 5 3d 1 4/ 1 


0.0711 


3s 1 3p 6 4s 1 


0.0849 


5 


3s 1 3p 5 3d 1 4p 1 


0.0639 


3s 1 3p 5 3d 1 4p 1 


0.0769 


6 


3s 2 3p 4 3d 1 4d 1 


0.0528 


3s 1 3p 5 3d 1 4f 1 


0.0738 


7 


3s 1 3p 5 4p 1 4d 1 


0.0517 


3s 1 3p 5 4s 1 4p 1 


0.0539 


8 


3s 2 3p 4 4/ 2 


0.0455 


3s 2 3p 3 3d 2 4p 1 


0.0535 


9 


3s 2 3p 4 4<2 2 


0.0399 


3s 2 3p 4 4/ 2 


0.0448 


10 


3s 2 3p 4 4s 1 4d 1 


0.0364 


3s 2 3p 4 3d 1 4d 1 


0.0441 


11 


3s 2 3p 5 4p 1 


0.0355 


3s 2 3p 4 4s 1 4d 1 


0.0413 


12 


3s 2 3p 4 4s 2 


0.0292 


3s 2 3p 4 4<2 2 


0.0375 


13 


3p 6 3d 2 


0.0272 


3p 6 3d 2 


0.0348 


14 


3s 1 3p 5 4d 1 4/ 1 


0.0261 


3s 2 3p 3 4p 3 


0.0298 


15 


3p 6 3d 1 M 1 


0.0232 


3s 2 3p 4 4s 2 


0.0291 


16 


3p 6 4p 2 


0.0207 


3s 1 3p 5 4p 1 4d 1 


0.0256 


17 


3s 2 3p 4 3d 1 4s 1 


0.0203 


3s 2 3p 3 3d 2 4f 1 


0.0252 


18 


3s 2 3p 4 4p 1 4/ 1 


0.0190 


3s 1 3p 4 4s 1 4p 2 


0.0234 


19 


3p 6 4s 2 


0.0178 


3s 1 3p 4 3d 2 4s 1 


0.0232 


20 


3p 6 4d 2 


0.0140 


3p 6 4p 2 


0.0209 


21 


3p 6 4/ 2 


0.0126 


3s 2 3p 4 4p 1 4/ 1 


0.0182 


22 


3s 1 3p 6 4s 1 


0.0046 


3p 6 4s 2 


0.0138 


23 






3s 1 3p 4 3d 1 4p 2 


0.0137 


24 






3p 6 4/ 2 


0.0122 


25 






3s 1 3p 4 3d 3 


0.0115 


26 






3p 6 3d 1 4d 1 


0.0107 


27 






3s 1 3p 4 3d 2 4d 1 


0.0087 


28 






3s 1 3p 4 4p 2 4d 1 


0.0070 


29 






3s 2 3p 4 3d 1 4s 1 


0.0069 


30 






3p 6 4d 2 


0.0044 


31 






3s 1 3p 5 4d 1 4f 1 


0.0026 


32 






3s 2 3p 3 4p 2 4/ 1 


0.0021 


33 
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Table 2. Values of p chosen to perform the corresponding open-core CI 
calculations for each MR, of CI and Cl — . 





MR 


selected p 


ci- 


MRD 
MRT 


1-6, 8, 10, 13, 14, 16, 19, 23 
1-9, 12, 14, 16, 20, 22, 27, 33 


Cl 


MR1 4 
MRI5 
MR2 


1-4 
1-5 
1-10, 12, 15 



and none in Is. These lists are finally reduced with respect to the full MR according 
according to the building rule (15 1. The resulting space is denoted MR-I/CVp. The p 
values selected for each MR ( 16||19 1 are presented in Table [2j 



4. Results and discussion 



4-1- Valence results 

The results of the MCHF and FC-MCHF calculations using the multireferences MRD 
for Cl~, and MR2 for CI, with n max l max = 9k, 10k and life are presented in Table [3] 
The 54 and Ssms values converge smoothly with the number of correlation layers 
but a significant difference is found between the MCHF and FC-MCHF approaches. 
The ground state of Cl~ is more stabilized by the optimization of the core orbitals 
than the one of CI. The impact on the difference of mass polarization parameters 
is enormous (« 36%), leading to a variation of the isotope shift of about 70%. The 
MCHF results of Table [3] are in agreement with the value of the differential mass 
polarization ASgMS = 0.091(25) <2q 2 extracted from the experiment [26] and the 
value pi). 

It should however be noted that, due to the total wave function symmetry 
properties, the core optimization in the MCHF calculations leads to significantly 
different contributions in Cl _ and CI. Indeed, the equivalence of the 2p and 3p 
subshells in the main configuration of Cl _ allows a large rotation of the P2 P {r) and 
Ps p (r) radial functions as it leaves the dominant contribution to the wave function 
(ls 2 2s 2 2p 6 3s 2 3p 6 ) invariant. Contrary to the anion, this effect is suppressed in the 
neutral atom by the smallness of the 2p e 3p 5 <H> 2p 5 3p 6 interaction. 

Figure [l] illustrates the dramatic mixing of core and valence occurring in Cl _ 
as compared to CI. The plot at the top shows the total radial density calculated 
with the program density [37] from the Cl~ MCHF wave function. In a closed core 
calculation, this total density can easily be separated into its [Ne] core contribution 
and its valence contribution. It is compared to the valence radial density of neutral 
chlorine. In the bottom panel of Figure [T] we present the total difference of radial 
density between Cl _ and CI as well as its core and valence contributions. It is usually 
accepted that the core could change only slightly when comparing the neutral atom 
with its negative ion pQ. This is indeed observed in the total radial density difference, 
the "extra" electron of Cl~ being efficiently screened out of the [Ne] core. However, 
we can also observe that it is the result of a large compensation between the core and 
valence contributions. It means that "freezing the core" does not mean the same thing 
in CI and Cl _ . This effect did not appear when comparing S~ to neutral sulfur [29 . 
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Table 3. Total energies (E, in Eh), Ssms parameters (in a 2 ) and size of the 
expansions (NCSF) in the MRD-I model of Cl~ and in the MR2-I model of CI for 
the largest active sets. The core-orbitals are either fixed to their HF shapes (FC- 
MCHF) or optimized (MCHF). The electron affinity (A) and A5 SMS calculated 
using these models are reported in the third part of the table. 









FC-MCHF 


MCHF 






nl max 


NCSF 


E 


Ssms 


E 


Ssms 


Cl-(MRD) 


9k 
10k 
Ilk 


113 691 
159 948 
214 416 


-459.838 5962 
-459.838 8245 
-459.838 9311 


-80.445 6256 
-80.445 8813 
-80.446 0807 


-459.839 4636 -t 
-459.839 6907 -i 
-459.839 7916 -i 


0.420 4925 
0.420 8569 
0.421 0037 


C1(MR2) 


9k 
lOfc 
llfc 


291 878 
410 462 
550 117 


-459.704 5818 
-459.704 7414 
-459.704 8192 


-80.324 8911 
-80.325 1537 
-80.325 3422 


-459.704 9422 -t 
-459.705 0994 -i 
-459.705 1767 -J 


0.331 9936 
0.332 2030 
0.332 3553 






max 


A 




A 




Diff 




9k 
Wk 
Ilk 


0.134 0144 
0.134 0831 
0.134 1119 


0.120 7345 
0.120 7276 
0.120 7384 


0.134 5215 
0.134 5913 
0.134 6149 


0.088 4989 
0.088 6538 
0.088 6483 


Exp. 










0.134 7(l) a 


0.091(25) 



a A™ (see Eq.[jj}. 



Figure 1. Top: total radial density of the Cl — MCHF wave function calculated 
with the program DENSITY 1371 and its [Ne] core and valence contributions. 
It is compared to the valence radial density of neutral chlorine. Bottom: 
total difference of radial density between Cl _ and CI and its core and valence 
contributions. 
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1 2 3 4 5 6 
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Table 4. Total energies (E, in Eh), $SMS parameters (in a 2 ) and size of the 
expansions (NCSF) in the MRT-I model for Cl~ , in the MRI4-I and MR1 5 models 
for CI. The active set is limited to 9k, 10k and life and the core orbitals are fixed 
to their HF shapes. 



FC-MCHF 





nl max 


NCSF 


E 


Ssms 




9k 


330 083 


-459.839 5099 


-80.432 1726 


Cr (MRT) 


10k 


467 617 


-459.839 7392 


-80.432 4083 




Ilk 


630 131 


-459.839 8463 


-80.432 6017 




9k 


65 286 


-459.703 5158 


-80.324 7147 


C1(MR1 4 ) 


10k 


91 915 


-459.703 6706 


-80.324 9490 




Ilk 


123 357 


-459.703 7494 


-80.325 1612 




9k 


152 483 


-459.704 4618 


-80.323 3889 


C1(MR1 5 ) 


10k 


212 132 


-459.704 6178 


-80.323 6181 




Ilk 


282 047 


-459.704 6972 


-80.323 8321 



The impact of the CI - main configuration symmetry on the optimal core orbitals 
is small enough that MCHF results can still be considered meaningful although a priori 
untrustworthy from a strict point of view. Nevertheless the lack of a clean distinction 
between core and valence contributions to the wave function becomes dramatic in an 
open-core CI approach. 

Table [4] presents the FC-MCHF calculations MRI4-I, MRI5-I of CI and MRT 
of Cl~ for [~9Af|, |"10fc] and [life]. Note that, from the comparison of Tables [3] 
and|3J the energy values of the MRI5-I and MR2-I models are in close agreement but 
the corresponding Ssms parameters do not compare so well, the mass polarization 
expectation values of the MR2 based approach being in better agreement with the 
MRI4-I model. It also should be observed that the MR2-I CSF spaces are almost 
twice larger than the ones of MRI5-I for a given number of correlation layers. 



4-2. Open-core CI analysis and error estimation 

In the following, we denote the CI and CI - models (r,p) and {r' ,p'), respectively, r 
denoting the used multireference, i.e. either l4,l 5 ,2 for CI or D,T for Cl~. The 
full array of results of the open-core CI calculations on CI and Cl~ are reported in 
Table [5j Considering that the results of the calculations performed on CI and Cl~ are 
independent would lead to separate studies of the convergence on each system and 
the summation of the so-deduced error estimations would give an overestimated final 
uncertainty on the theoretical e A and corresponding IS. This procedure indeed does not 
take account of the fact that the correlation contributions to the separate systems tend 
to cancel out when estimating the differential effects. To extract the differential ASsms 
value from our results, we need some guideline. In valence calculations, this guideline 
is the number of correlation layers. Unfortunately, in the open-core CI approach, there 
is no satisfactory intrinsic convergence path in the space of the calculation parameters 
r, r', p and p' . Thanks to the correlation observed between the energy and the 
mass polarization expectation value [2H US] , we use the reference e Af^ estimated in 



Section 2.1 as a guideline. 
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Table 5. Full array of open-core CI (r,p) results for CI and CI (see text). All 
values are given in in atomic units. 
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-459.7686643 


-79.756779 


3 


-459.9015022 


-79.876465 




5 


-459.7688295 


-79.758733 


4 


-459.9018425 


-79.880946 




6 


-459.7691030 


-79.758821 


5 


-459.9024779 


-79.870426 




7 


-459.7692978 


-79.757749 


6 


-459.9028828 


-79.869798 




8 


-459.7693594 


-79.757970 


7 


-459.9031145 


-79.870329 




9 


-459.7694374 


-79.756594 


8 


-459.9032110 


-79.871590 




12 


-459.7695755 


-79.752764 


9 


-459.9035007 


-79.863268 




15 


-459.7695952 


-79.752384 


12 


-459.9038527 


-79.857397 










14 


-459.9039916 


-79.855950 










16 


-459.9040884 


-79.854589 










20 


-459.9042568 


-79.850471 










22 


-459.9042933 


-79.850431 










27 


-459.9043424 


-79.849504 










33 


-459.9043545 


-79.849094 
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Figure 2. E versus Ssms ( m a.u.) plots of the results of the open-core CI 
calculations (see Table [5| for CI - (a) and CI (b). Convergence is read from right 
to left (decreasing energy) . The plain line of Figure (a) is the relation \21\ and 
the dashed lines correspond to the uncertainty on the fitted parameters (sec text). 




-459.905 -459.902 -459.899 -459.771 -459.768 -459.765 

E(CT) (E h ) E(CI) (E h ) 



Figure [2] plots the Ssms results obtained in the (r,p) model (5 , sMs( r iP)) I versus 
the energies E(r,p), for CI and CI - in Figures (a) and (b) respectively. The p- values 
are given in Table [2] Note that the same scale is used in both plots. The p = 1 
points are outside the frames. In order to extrapolate the results of the MRI5 model, 
we add the configuration 3s 1 3p 4 4s 1 4p 1 to MR1 5 , generate the corresponding open- 
core MRlg-I/CVg [life] CSF list and perform a CI calculation on this space using the 
MRI5-I orbital set. 

Figure [2]-a exhibits a good agreement of the results using the two references of 
Cl _ . Furthermore, the convergence inp (from right to left in the plots) is remarkably 
clean. Even though the correlation between E(r,p) and S'sms ( r 'P) °f ^1 calculations 
is apparent in Figure [2}b, it is less precise. 

In these two plots, we see that (H 2 /m e ) S'sms is much more sensitive to p than the 
total energy. Furthermore, we see from the comparison between Tables [3] and |4j and 
Table [5] that the Ssms parameter is much more affected than the total energy E when 
open-core configurations are included in the expansion. This observation is supported 
by the larger impact of the core relaxation on the Ssms than on the energy in valence 
MR-I models (see Table [3]) . Fitting the eight last points based on the MRT model of 
Cl~ (see Tables [5]) , one obtains (in atomic units) 

SsMs(Cr) = -16.6(7) [E(Cr) - E] - 79.8538(2) (21) 

where E = —459.904 085 is the average of the energies of the fitted points. If all 
the aligned points of Figure [2]-a are fitted, i.e. if the five last points for MRD are also 
included in the fit, the angular coefficient becomes —17.5(4). It gives an indication 
of the uncertainty of the fit considering the convergence with respect to the used MR 
(r'). The relation (21) is illustrated by a plain line in Figure [2ja, the dashed lines 



being the lower and upper limits for an 1.3 uncertainty o n th e angular coefficient 
and a 2.10~ 4 <Zg 2 uncertainty on the offset. Using equation ( pi] ) with these upgraded 
uncertainties allows to extrapolate the Cl~ energy and corresponding Ssms 

ASsms = Ssms (CI) + 16.6(13) [E(C\) - E - e A^f] + 79.8538(2) 

(22) 

taking e AfJ^ = 0.1347(1) E/j from eq. The convergence of the total weight of the 
p first configurations of MRT-based wave functions with the corresponding energies 
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Table 6. Total energies and SgMS parameters of CI and ASgMS values 
extrapolated from the (r,p) results of CI to have E(Cl~) = E(Cl) - e A^Jj . 
•SsmsCCI - ) and ASgMS are deduced from equation \22\ . All values are in atomic 
units. 



CI model CI extrapolation 



T 


I) 




So M c:(Cl _ )far 2 l 


f\ tS'anAa \n.~ ^1 
*— ""oMd l J 


U 


1 


-459.8985 


-79.9456(98) 


0.1038(98) 


u 


1 


-459.8994 


-79.9312(86) 


0.0897(86) 


u 


2 


-459.9013 


-79.8995(59) 


0.1302(59) 




3 


-459.9020 


-79.8884(50) 


0.1201(50) 


Is 


2 


-459.9022 


-79.8840(46) 


0.1203(46) 


U 


4 


-459.9026 


-79.8785(41) 


0.1179(41) 


U 


3 


-459.9030 


-79.8719(36) 


0.1096(36) 




4 


-459.9036 


-79.8619(27) 


0.1074(27) 




5 


-459.9039 


-79.8567(23) 


0.1019(23) 


is 


6 


-459.9044 


-79.8478(25) 


0.0955(25) 


2 


1 


-459.8986 


-79.9437(96) 


0.1237(96) 




2 


-459.9015 


-79.8960(56) 


0.1602(56) 




3 


-459.9028 


-79.8753(39) 


0.1132(39) 




1 


-459.9034 


-79.8641(29) 


0.1073(29) 




5 


-459.9036 


-79.8614(27) 


0.1026(27) 




6 


-459.9039 


-79.8568(23) 


0.0980(23) 




7 


-459.9041 


-79.8536(21) 


0.0958(21) 




8 


-459.9041 


-79.8526(22) 


0.0946(22) 




9 


-459.9042 


-79.8513(23) 


0.0947(23) 




12 


-459.9043 


-79.8490(24) 


0.0962(24) 




15 


-459.9044 


-79.8487(25) 


0.0963(25) 



suggests that, down to an energy of —459.9052 E/,., the variational principle is not 
violated within the limits of our model. 

Table [6] gives the extrapolated energies and Ssms °f ^ ne d an< i ^ ne resulting 
^•S'sms values corresponding to each open-core CI calculation performed on CI (see 
Table [5]). We regroup the results of the models with r = I4, 15 and 1§ and order 
them by decreasing energy. The quoted errors only take into account the uncertainty 
deduced from the Cl _ calculations convergence and from the estimated uncertainty 
on 

In the r = 2 calculations, the increase of the Ssms (CI) with p is close to the 
angular coefficient of equation (21 1 giving Ssms (CI - ). It leads to the observed 
smoother convergence of the r — 2 results of Table [6] than the ones based on the 
other combined r = 14,15, lg approaches. Still, all the calculations agree with the 
experiment (AS SM S = 0.091(25) a^ 2 ). From Table^ we adopt the value 0.096(3) a^ 2 . 
Using the fact that the results of (r.p) — (2,15) and (16,6) models are consistent, 
we estimate the uncertainty associated to the convergence of the calculations on CI 
by taking the difference between the (2,12) and (ls,5) values, leading to the final 
estimation ASsms = 0.096(9) a^ 2 . 
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Table 7. Experimental and theoretical isotope shifts (IS) on the e A for the 
chlorine isotopes A = 37,35. The experimental NMS of 0.741 GHz is used to 
obtain the IS. All values are in GHz. Valence MCHF and FC-MCHF are based 
on the models MRD and MR2 respectively for Cl — and CI. The final value is 
deduced as explained in Section [4^2] 



SMS 



MS 



FS 



RIS 



IS 



HF 

val. FC-MCHF 
val. MCHF 
final results 



Exp. 
DF 

MB low corr. 



-1.348 
-0.674 
-0.495 
-0.535(51) 



-1.3 
-0.50 



-0.607 
+0.067 
+0.246 
+0.206(51) 



-0.6 
+1.24 



This work 

-0.003(22) 
-0.002(20) 
-0.003(21) 
-0.003(22) 

Berzinsh et al |26| 



-1.351(22) 
-0.676(20) 
-0.497(21) 
-0.538(72) 



1-0.014(14) 
1-0.014(14) 



-0.51(14) 

-1.3 
+0.51(2) 



-0.610(22) 
+0.065(20) 
+0.244(21) 
+0.203(72) 

+0.22(14) 

-0.6 
+1.26(2) 



5. Comparison with experiment and conclusion 



In Section |33| we elaborate large scale valence MCHF models and subsequently add 
open-core CSFs expansions. Focusing on differential effects, more flexibility is added 
in the CI - models than in the ones of CI. The correlation between the energy and 
the mass polarization expectation value convergences in a sequence of calculations, 
already observed in earlier work [49], allows to do realistic extrapolations around the 



results of the largest Cl~ calculations (see Section 4.2 1. We hence fit our results to 
the experimental non-relativistic electron affinity estimated in Section [2~Tj A robust 
uncertainty for the (valence+core-valence) correlation contributions to the isotope 
shift is deduced. However, core-core correlation effects are known to be important 
in isotope shifts calculations [50] EI] and more accurate approaches should be used, 
ultimately including relativistic corrections. 

Table [7] compares results of our work with the experimental and theoretical values 
of Berzinsh et al |26j?) The experimental normal mass shift (NMS) of 741 MHz is used 
to deduce the total IS so that a cleaner theory/experiment comparison is found in the 
residual isotope shift (RIS = IS — NMS). The uncorrelated Dirac-Fock (DF) results 
from Berzinsh et al [26] are also reported. For each model, we give the specific mass 
shift (SMS), total mass shift (MS = NMS + SMS) and field shift (FS). As pointed out 
by Berzinsh et al [55] , their estimated field isotope shift of the electron affinity is well 
below the error bars of the IS measurement. The difference between their FS values 
and ours is caused by another choice of 8 (r 2 ) change between isotopes (see discussion 



in Section 2.3). With this respect, the uncertainty associated with the very small FS 
constitutes a large source of uncertainty in the final RIS and IS values. 

The calculations of Ref. [55] include the relativistic effects but only the lowest 
order correlation corrections, predicting a RIS opposite to the measured one. It is 
clearly seen from the FlG. 6 of this article displaying the behaviour of the SMS for 
different levels of approximation, that the theory-observation discrepancy is due to an 
insufficient treatment of the correlation effects. 

§ Note that in their equation (10) expressing the separation between the mass shift and the field 
shift, Berzinsh et al adopted the first term of jll| for the mass contribution. 
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4.1 



The results of the MRD-I|~llfc] for CI" and MR2-I|~llfc] obtained in Section 
are combined to give the valence frozen core (FC-MCHF) and optimized core (MCHF) 
estimations presented in Table [7] At first sight, relaxing the core orbitals lifts the 
discrepancy between the valence model and the experimental values. However, in 
this work, we observe that the optimization of the core orbitals of CI and CI - leads 
to significantly different contributions, a priori destroying the balance between the 
core/valence distinction of the two systems. This effect arises from the occupation 
numbers symmetry in the CI - main configuration. If it might be a good reason to 
distrust the closed-core MCHF values of Table [7] it is not a demonstration that they 
should be ignored. 

The isotope shift on the electron affinity is found to be anomalous for 
beryllium [52] and oxygen |16[ 128] , corresponding to a smaller electron affinity for the 
heavier isotope. For sulfur, the theoretical estimations also predict an anomalous IS, 
allowed by the large experimental error bars, although the latter are centered on 
a positive difference ( 54( 34 S) - e A( 32 S) = +0.0023(70) cm" 1 ) [25]. The IS on the 
detachment thresholds were also found to be negative in carbon [3D], waiting for 
some experimental confirmation. The isotope shift on the chlorine electron affinity is 
definitely normal, the NMS largely dominating the SMS. For all these systems, the 
FS are at least two orders of magnitude smaller than the total mass shift, to put 
in contrast with heavy negative ions, such as Os~ for which the field shift clearly 
dominates the IS on the electric-dipole bound-bound transition frequency [53]. For 
systems for which the FS can be neglected, the isotope shift on the electron affinity 
can be rationalized by looking at equation (11). The electron affinity being positive 
for a bound negative ion, the isotope shift will be normal if 54 > (h 2 /m e )ASsMS- 
Taking the atoms along the third period, from Al to CI, as an example, the electron 
affinity increases regularly from the left to the right of the periodic table, with a dip 
corresponding to phosphorous, while the ASsms positive values defined by eq. Q are 
relatively stable. Aluminum and silicon having a rather small 54, they are expected to 
have a significant negative (anomalous) isotope shift while chlorine, as demonstrated 
in the present work, is characterized by a large electron affinity that determines the 
(normal) IS. The normal-anomalous IS change corresponds to the crossing between 
the (ti 2 /m e )ASsMS an d 54 trends occuring around sulfur (Z = 16) which IS sign on 
the electron affinity is not fully definite [29] . 

As pointed out by Andersson et al [20J, the analysis of the electron affinities 
of atoms shows that, with the exception of the hydrogen atom, calculated electron 
affinities are less accurately known than measured ones. However, we have shown 
in [2D] that the theoretical IS values on the electron affinity could be more accurate 
than the observed ones. In the present work, we demonstrated that theoretical 
calculations of the electronic factors of the mass contribution have evolved remarkably, 
solving a longstanding discrepancy between calculated and measured specific mass 
shift on the electron affinity of chlorine. 
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